Peierls-type Instability and Tunable Band Gap in Functionalized Graphene 
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Functionalizing graphene was recently shown to have a dramatic effect on the electronic properties 
of this material. Here we investigate spatial ordering of adatoms driven by the RKKY-type inter- 
actions. In the ordered state, which arises via a Peierls-instability-type mechanism, the adatoms 
reside mainly on one of the two graphene sublattices. Bragg scattering of electron waves induced 
by sublattice symmetry breaking results in a band gap opening, whereby Dirac fermions acquire a 
finite mass. The band gap is found to be immune to the adatoms' positional disorder, with only an 
exponentially small number of localized states residing in the gap. The gapped state is stabilized 
in a wide range of electron doping. Our findings show that controlled adsorption of adatoms or 
molecules provides a route to engineering a tunable band gap in graphene. 



The unique electronic properties of graphene, a one- 
atom-thin carbon sheet with a tunable electron den- 
sity in and high carrier mobility 0, 0], make it an at- 
tractive material for applications in nano-electronics 
However, because of the gapless semi-metallic charac- 
ter of graphene band structure, the future of graphene 
electronics depends on developing methods to engineer a 
band gap in this material. The gapless character of elec- 
tron dispersion in pristine graphene is protected by the 
high symmetry of its lattice, in which two carbon sites 
in the unit cell are equivalent. The simplest kind of gap- 
opening perturbation which lifts this symmetry can be 
described by unequal potentials ua and Ms on the A and 
B sites [^, leading to a finite mass of Dirac quasiparti- 
cles near points K and K' of the Brillouin zone. The 
quasiparticle spectrum, described by the Hamiltonian 



H 



K(K') 



vo « lO'^m/s, 
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features a band gap of size A = \ua — ub\, which opens 
due to Bragg scattering of electron waves on the periodic 
sublattice potential. 

A gap opening via such a mechanism could occur in 
epitaxial graphene, grown or placed on a lattice-matched 
substrate |,6|~i8i]. Yet, while the approach involving lattice- 
matched substrates is simple and direct, combining it 
with transport measurements proved challenging (see 
also Ref. [9|). A gap opening due to sublattice asymme- 
try is more readily achievable in bilayer graphene, where 
the sites A and B reside on different layers. In bilayer 
graphene, the A/B symmetry can be lifted by asymmet- 
ric chemical doping or electrical gating 12 1, lead- 
ing to a gap opening. 

Another promising method for gap engineering relies 
on spatial confinement, involving patterning graphene 



into narrow ribbons |13l . |14| , or quantum dots |15| . The 
gap obtained by such a method can be tuned by varying 
spatial width of graphene ribbons or dots. However, the 
approaches relying on spatial confinement are prone to 
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FIG. 1: Peierls-type instability resulting from adatom or- 
dering over sublattices A and B. A gap in the density of 
electronic states opens up due to Bragg scattering on the 
A/B modulation when the occupation probabilities are un- 
equal, Ua 7^ riB. The states in the gap move down in energy 
into a peak positioned at the energy of a single adatom reso- 
nance. Inset: the ordered state is stabilized in a wide range 
of carrier densities, for which the energy gain per adatom is 
positive. For details of calculation see discussion following 
Eq.®. Different curves correspond to the occupancy frac- 
tion UA/iuA + ub) = 0,0.1, ...,0.5 with {ua + nB)/2 = 0.04, 
U = 6to, W = ito. 

disorder, because of scattering of electron waves on rough 
edges of patterned graphene. Localized states, appearing 
inside the band gap, transforms it into a "transport gap" 
[l6| . In contrast, the gap opened due to lifting the A/B 
symmetry can be expected to be more robust in the pres- 
ence of disorder, as long as the mean free path is large 
compared to the A/B modulation period. 

An elegant approach to modify electronic properties of 
graphene, demonstrated recently 17, l^j, is based on the 
well established technique of chemical functionalization. 
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in which groups such as H, OH, or F bind covalently to 
carbon atoms, transforming the trigonal sp^ orbital to 
the tetragonal sp^ orbital. Such transformation drasti- 
cally alters local electronic properties. Theory predicts 
that, at 100% coverage by H adatoms, graphene turns 
into a wide-gap semiconductor, called graphane [l^. The 
experiments however, are done at low coverage, 

typically of about few percent. Can a state with a band 
gap be realized in the low-coverage regime? 

Electronic properties at low adatom coverage are dom- 
inated by resonant scattering of electron waves on the 
adatoms Pairwise RKKY-type interactions be- 

tween adatoms were analyzed in Ref. [23| . where the in- 
teraction sign was found to depend on whether the in- 
teracting atoms occupy the same sublattice or different 
sublattices. Such sublattice dependence suggests that the 
RKKY interactions can drive ordering of the adatoms in 
which sublattices A and B become unequally populated. 

Here we propose a mechanism for spontaneous order- 
ing, illustrated in FiglU which is analogous to that of 
Peierls instability. The adatom ordering over sublattices 
A and B leads to a gap opening due to electron waves 
Bragg scattering on the A/B modulation, resulting in 
electronic states in the gap shifting up and down in en- 
ergy. Crucially, these shifts are asymmetric, with states 
shifting predominatly down in energy to a peak centered 
at the energy of a single atom resonance, e = £o < 0. 
The system gains energy as a result of such level shifts 
for electron dopings in the range indicated in Fig [T] inset, 
corresponding to positive chemical potential values. For 
such dopings, the gapped state with unequal sublattice 
population is stabilized. The gap value is determined by 
the scattering properties of adatoms and their concentra- 
tion, and is therefore tunable. 

Because of the resonant character of electron scatter- 
ing, the electron-mediated interactions fall off slowly with 
distance at adatom separations r ^ Iq = hvo/lsol, as 
U{r) ~ 1/r, and more rapidly at larger distances [23|. 
Hence for not too low adatom coverage, n ^ n* = 
(a/£o)^j where a = 0.142 nm is the lattice constant, the 
adatom ordering cannot be analyzed using a pairwise in- 
teraction model. Here we present a theory which fully ac- 
counts for the non-pairwise, collective nature of electron- 
mediated interactions in functionalized graphene. 

Our approach applies to different atoms and chem- 
ical groups used to functionalize graphene. Ab initio 
study |22l | predicts the resonance energy values which 
span a wide range: —Eq = 0.03, 0.11, 0.70, 0.67 eV for 
H, CH3, OH, and F, respectively. This corresponds to 
the characteristic values ss 10""*, 10"'^, 0.05, 0.05. In 
the limit of very small coverage, n ^ n*, a pairwise in- 
teraction model can be used to describe ordering 24 , 2^ , 
whereas for larger coverage values a self-consistent treat- 
ment presented below must be employed. 

The Peierls-type scenario for ordering described above 
can be tested by direct numerical diagonalization of the 
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FIG. 2: The density of electronic states as a function of en- 
ergy for different sublattice occupancy ratios, obtained by nu- 
merical diagonalization of the Hamiltonian ((2} , averaged over 
160 realizations of disorder. The band gap, which opens at 
7^A ^ nB and ub <^ ua, is immune to disorder: no electronic 
states are found inside the gap. Sublattice ordering results in 
the energy gain at positive dopings (inset). Parameters used: 
system size 62 x 34, [/ = 6fo, n = [ua +nB)/2^ 0.034. 

nearest neighbor tight binding Hamiltonian (see Figl2|): 



to(V'i^x'+h.c.)-f ^M(x)7/.t^, 



(2) 



-x'| = l 



with to ~ S.leV, and potential u(x) ~ J2iUS{:x. — x^) 
taking value U on the sites occupied by adatoms. Large 
U :s> to was used to model the effect of the sp^ to 
sp^ transformation, which inhibits the conduction elec- 
trons from occupying the adatom sites, effectively turn- 
ing these sites into vacancies. In the simulation shown in 
Figl2]we used /7 = 6to, which gives the resonance energy 
positioned at £0 ~ — 0.4eV. 

The behavior of the DOS, obtained for different oc- 
cupancy values by diagonalizing the Hamiltonian ([2]), 
agrees well with the results obtained by an analytic 
method (see Fig[T]). The peak at e = SQ) which is a signa- 
ture of resonant scattering on individual adatoms (20l . ]2^ , 
is present for all occupancies, but is more pronounced for 
the sublattice-ordered state, ^ or n^. The 

resonances marked 2' and 2" correspond to the single- 
particle states formed near two neighboring adatoms. 
The DOS remains finite at all energies for ua ~ ub- In 
contrast, the DOS vanishes in the interval < e < QAeV 
for the sublattice-ordered state, which corresponds to the 
band gap opening. 

To estimate the energy gain due to ordering, we eval- 
uate the energy of the system as a function of the 
adatom occupancy fraction and electron concentration 
N/Ns — 1/2 {N is the number of electrons, Ng is the 
number of sites). The results for 5E^ the energy gain 
per adatom, shown in Figl2]inset, agree with our Peierls- 
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FIG. 3: Two phases of adatoms on graphene: disordered (a), 
where adatoms are randomly distributed over two sublattices 
of the hexagonal lattice, A and B, and ordered (b), where 
adatoms preferentially occupy one of the sublattices. The 
adatoms residing on sublattice A{B) are shown in red (green). 
Schematics of the energy spectrum in the two phases are also 
shown. In the ordered phase, the sublattice symmetry break- 
ing leads to a band gap. 

type argument: the sublattice-ordered state is stabilized 
at positive doping. The ordering temperature, estimated 
as Tc — 2SE (see Appendix), takes values in the hundreds 
of Kelvin for parameters used in Fig. [21 

Remarkably, in our simulation the energy gap is found 
to be immune to the effects of disorder. As Figl2] in- 
dicates, no states inside the gap are found when the 
adatoms randomly populate one sublattice, A or B (160 
disorder realizations were analyzed). Such partial order- 
ing, illustrated in Fig|31 is sufficient to completely expel 
the midgap states, and open a band gap. This surprising 
result seems to be in an apparent contradiction with the 
intuition based on the physics of localization of electronic 
states in disordered systems. 

To understand this behavior, we shall start with a 
simple case of a weak adatom potential U <^ to, and 
then generalize to the case of strong potential. At a 
weak potential, within the mean-field approximation, we 
have ua = uaU and ub = ubU in the Hamiltonian 
(HI). For the gap to survive in the presence of disor- 
der, its value should exceed the disorder broadening, 
A = \ua — ub\ ^ h/r^r^A- Using the Born approximation 
for the scattering rate, we find 

no = „ ^ (3) 



-J^ = :^pie)\U\^inA+nB), 
T(e) 2no 



where 2no is the density of carbon atoms in the graphene 
lattice. Taking p{e) to be that of pure graphene for one 
spin projection, p[e) = ^^2^2 ; and using the relation vq — 
Moa/2h, wc rewrite the condition A 3> ft/re^A as 

{nA + nB)U^ <t:tl. (4) 

This condition is always satisfied for weak adatom po- 
tential U ^ to. 



This argument can also be applied, with a slight mod- 
ification, to the case of strong adatom potential, U to- 
This can be done by replacing U in Eq.Q by a suit- 
ably defined T-matrix, see Eq.([5]). This again yields the 
condition Q, which is satisfied at low enough adatom 
coverage. Furthermore, as numerical results presented 
in Fig. [2] show, the gap persists even at higher adatom 
concentration, [tia + nB)U^ ~ to- 

Of course, strictly speaking, the DOS inside the gap 
must be nonzero. However, since the states deep in the 
gap can arise only due to relatively large fluctuations of 
disorder, their contribution to the DOS is exponentially 
small. To analyze this quantitatively, we shall focus on 
the simplest case of a weak potential, J7 <C toi assuming 
that all adatoms reside on the A sublattice. We model the 
effect of disorder by the Hamiltonian ([ij with fiuctuating 
gap A(r) = UA = nA{'r)U, and zero ub- 

This problem can be mapped on the well studied prob- 
lem of the DOS below the band edge in a disordered 
semiconductor [2^, [23|. Starting with the equations 
eipA — A(r)^/;^ +voP+'ipB, £'ipB — voP-'^pA, and eliminat- 
ing the component 'lpB^ we obtain an eigenvalue equation 



(eA(r)-t-2v2)^^(r) = £VA(r). 



(5) 



For e near the upper band edge, given by the disorder- 
averaged potential, A = uaU , we expand in (5e = e — A 
to bring the eigenvalue equation to the form of the 
Schroedinger equation for a massive non-relativistic par- 
ticle: 



(5£-0A(r) 



V2 



(6) 







where (5m (r) = 5n^(r)C/ is the fluctuating part of the 
gap A(r). We treat the long wavelength fluctuations of 
A(r) as Gaussian with the two-point correlation func- 
tion {Su{r)Su{r')) = U^{6nAir)SnAir')) = XS{r - r'), 
where A = C/^n^i/no, as appropriate for delta-correlated 
adatom positions. In this case the DOS for the problem 
(|6]) decays exponentially away from the band edge [27| . 

p(0 < £ < Ua) oc exp(— c|fc|/i^), ^ = Xm (7) 

with c a constant of order one. Estimating the energy 
scale we see that, at low coverage ^ 1, it is much 



smaller than the gap width: ^ = Xm 



2novfi to 

A. Therefore we conclude that the DOS is exponentially 
small within the gap. This is consistent with the results 
of our simulation, Figl^J in which some smearing of the 
DOS was observed at the gap edge; however no states 
were found deep inside the gap for all of the 160 disorder 
realizations which we analyzed. 

An analytic approach to analyze transport properties, 
such as disorder scattering and conductivity, can be de- 
veloped using a self-consistent T-matrix approximation 
(SCTA) (see Appendix). In the model 1^, the T-matrix 



4 



of an individual adatom has a resonant form , 



U^U/no, (8) 



e\n{iW/e) + S' U 

where W ~ 3^0 is the bandwidth, and the parameter 5 
determines the energy of the resonance, eolnW/jeol = 
~5. We note that the description of adatoms by an on- 
site potential is equivalent to the model which describes 
adatom states in terms of a localized level hybridized 
with the graphene continuum 22| , since the form of the 
T-matrix at low energies is the same in both approaches. 

We calculate disorder-averaged Greens functions using 
the self-consistent approach (see Appendix). The Greens 
functions are then used to extract the density of states, 
p{e) = — ilmTrG'(e + iQ). The resulting energy depen- 
dence of the DOS is shown in FigHJ The main features, 
such as the gap opening for imbalanced A/ B occupan- 
cies, accompanied by the shift of electronic states into the 
resonance peak of a single adatom, are in agreement with 
our numerical results. The DOS is then used to evaluate 
the total energy as a function of the chemical potential, 
= ep{e)de. The change in the electronic energy 
due to sublattice ordering, 5E = i?„^=„B — Enj^^uB^ 
depends on /i as shown in FigH] inset. The sublattice- 
ordered state is stabilized for positive dopings > 0. 

The gap opening and its character reveals itself in 
transport measurements, since the temperature depen- 
dence of conductivity is activation-like in systems with 
an intrinsic band gap, a cx ^-^/f^^T ^ j^^^g variable- 
range hopping behavior of the Mott or Efros-Shklovskii 
form for systems with a transport gap. In addition, sub- 
lattice ordering suppresses scattering, which leads to an 
increase in conductivity for electron doping above or be- 
low the gap 28^. Observation of such an increase can 
serve as a hallmark of adatom ordering. 

Conductivity of the system, evaluated within 
SCTA approach, is given by 28|, 

' ■ m ■ 



the 



-K^h? vl{hA + ub) n'^h {ha + ns)^ 1^1 

» A, 



(9) 



l^ol 



which is valid far outside the gap region, 
The quadratic dependence on p, in the leading term is 
characteristic for resonant scattering js^l- The sublead- 
ing term, proportional to UAns / (riA +nB)'^, is negative. 
Hence, scattering is indeed suppressed and conductivity 
enhanced due to sublattice ordering. 

An experimental approach to realizing the sublattice 
ordered state depends on the lateral mobility of adatoms. 
If the adatoms remain mobile below the ordering temper- 
ature Tc — 26 E (see Appendix), the ordering will occur 
via a conventional Ising-type phase transition. In this 
case, a rapid cooldown following ordering may be needed 
to prevent adatom clustering [23|. However, if the mobil- 
ity is quenched at temperatures T > Tc, the system must 
be annealed at T ^ Tc to achieve ordering and gap open- 
ing. Since only a small part of adatom's entropy needs to 



be removed for sublattice ordering (see FiglS]), it should 
take only a few hops by each adatom to transition into 
the gapped state. 

In summary, the interaction between adatoms in func- 
tionalized graphene can drive sublattice ordering via a 
Peierls-type transition. The band gap, opened by Bragg 
scattering of electron waves on the sublattice modula- 
tion, is immune to positional disorder of adatoms, with 
the density of localized states inside the gap being expo- 
nentially small. The gapped state is shown to be stable 
in a wide range of electron doping. 

We thank A. K. Geim, K. S. Novoselov, and E. Roten- 
berg for useful discussions. This work was supported by 
Office of Naval Research Grant No. N00014-09-1-0724. 



[1] 
[2] 

[3] 
[4] 

[6: 

[7 

[s: 

[9 

[lo; 

[11 

[12 
[13 

[14 

[15 
[16 

[ir 

[18 
[19 

[20 
[21 

[22 

[23 

[24 

[25; 

[26 
[27 
[28' 
[29 
[30 



K. S. Novoselov et al, Science 306, 666 (2004). 

J. H. Chen et al, Nature Nanotech. 3, 206 (2008). 

K. I. Bolotin et al, Solid State Comm. 146, 351 (2008). 

A. K. Geim, K. S. Novoselov, Nature Mater. 6, 183-191 

(2007). 

A. H. Castro Neto et al. Rev. Mod. Phys. 81, 109 (2009). 

C. Berger et al, J. Phys. Chem. B 108, 19912 (2004). 
S. Y. Zhou et al, Nature Mater. 6, 770 (2007). 

G. Giovannetti, P. A. Khomyakov, G. Brocks, P. J. Kelly, 
and J. van den Brink, Phys. Rev. B 76, 073103 (2007). 
E. Rotenberg et al. Nature Mater. 7, 258 (2008); S. Y. 
Zhou et al, ibid. 259 (2008). 

T. Ohta, A. Bostwick, T. Seyller, K. Horn, E. Rotenberg, 
Science 313, 951 (2006). 

E. V. Castro et al, Phys. Rev. Lett. 99, 216802 (2007). 
J. B. Oostinga et al. Nature Mater. 7, 151 (2008). 
Y.-W. Son, M. L. Cohen, and S. G. Louie, Phys. 
Lett. 97, 216803 (2006). 

M. Y. Han, B. Ozyilmaz, Y. Zhang, P. Kim, Phys. 
Lett. 98, 206805 (2007). 

L. A. Ponomarenko et al, Science 320, 356 (2008) 
S. Adam, S. Cho, M. S. Fuhrer, and S. Das Sarma, 
Rev. Lett. 101, 046404 (2008). 

D. C. Elias et al. Science 323, 610 (2009). 

A. Bostwick et aZ.. larXiv: 0904.2249 , unpublished. 

J. O. Sofo, A. S. Chaudhari, G. D. Barber, Phys. Rev. B 

75, 153401 (2007). 

V. M. Pereira et al, Phys. Rev. Lett. 96, 036801 (2006). 
V. M. Pereira, J. M. B. Lopes dos Santos, A. H. Castro 
Neto, Phys. Rev. B 77, 115109 (2008). 
T. O. Wehling, M. I. Katsnelson, A. I. Lichtenstein, Phys. 
Rev. B80, 085428 (2009). 

A. V. Shytov, D. A. Abanin, L. S. Levitov, Phys. Rev. 
Lett. 103', 016806 (2009). 

V. V. Cheianov, V. I. Fal'ko, O. Syljuasen, and B. L. 
Altshuler, Sol. State Comm. 149, 1499 (2009). 
V. V. Cheianov, O. Syljuasen, B. L. Altshuler, and V. I. 
Fal'ko, Europhys. Lett. 89, 56003 (2010). 

B. Halperin, M. Lax, Phys. Rev. 148, 722 (1966). 
Zittartz, J. S. Langer, Phys. Rev. 148, 741 (1966). 
A. Abanin, L. S. Levitov, arXiv: 1008. 1424 (2010). 
Pepin, P. A. Lee, Phys. Rev. B 63, 054502 (2001). 
M. Ostrovsky, I. V. Gornyi, and A. D. Mirhn, Phys. 



Rev. 



Rev. 



Phys. 



Rev. B 74, 235443 (2006) 



5 



APPENDIX 
Self-consistent T-matrix approximation 

Here we discuss the self-consistent T-matrix approxi- 
mation (SCTA). We use the Greens function expressed 
through disorder-averaged self-energy, 



G(e,p) 



-ip £b 



(10) 

with = <o(l + e ^^^^+e *''''^). An infinitesimal imagi- 
nary part ±iO should be added to e to obtain the retarded 
and advanced Greens functions. 

The self-energy is approximated by the average value of 
the T-matrix, separately for the A and _B-type adatoms, 

SA(e) = nA{TA{£)), ^b{£) = ?is(rs(e)), (11) 

with nji(^B) = 'T'A(B)'^o the adatom densities. Here the 
quantities 7^(5) i written as a 2 x 2 matrix, are given by 



Ta 
Tb 



U 



1-Ug 



9 = 



G(e,p). (12) 



The integral of the Greens function over the Brillouin 
zone is dominated by the regions near K and K\ giving 



9 



H-W^/saSb) 



SB 
EA 



W w 3to- (13) 



valid for |eyi(B)| <C W. Combining this result with 
Eg. dlip . we obtain two coupled equations for e^, eb'- 



EA 



uaU 



Eb 



ubU 



1 + IEb ' " 1 + IE A ' 

where we used the relation U = Uuq, and defined 
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1 U 



In 



eaEbJ 



(14) 



(15) 



To obtain retarded (advanced) Greens function, one 
should choose the branch of the logarithm that is ana- 
lytic in the upper (lower) half-plane. Solving numerically 
for EA, Eb as a function of e, we find the Greens function 
(fTO| and use it to calculate the density of states. 



p{e) = ilmTrG(e-l-iO). 

TT 



(16) 



This approach was used to produce the curves shown in 
Fig.l of the main text. 

An estimate of the ordering temperature 

In this section, we consider adatom ordering on the 
sublattices A and B, and estimate the ordering temper- 
ature. Since the ordering involves breaking of the sub- 
lattice symmetry, it can be described as a second-order 



phase transition. The broken symmetry is Z2, corre- 
sponding to a transition of the conventional Ising type. 

We define the order parameter in terms of ua and Ub, 
the populations of the A and B sites, as the imbalance 
in sublattice population: 



Sn= ^((n^)- {tib}), 



(17) 



where (...) denotes ensemble average. We expect that the 
order parameter Sn is zero above the transition tempera- 
ture, T > Tc, becoming nonzero at T < Tc- The behavior 
in the vicinity of the transition is governed by the depen- 
dence of the free energy F = E — TS. Below we analyze 
the dependence of F vs. Sn, and use it to estimate the 
ordering temperature Tc- 

The free energy of a generic state of adatoms is given 
by a sum of the interaction energy, analyzed in the main 
text, and the entropic contribution due to adatoms pop- 
ulating the A and B sites, 

N N 
FinA^ns) = Eint{nA,nB) - —Ts{nA) - —Ts{nB), 

(18) 

where N is the total number of carbon atoms, and 
s{nA(B)) is the entropy per site, which is given by 



s{x) = ~{xhix + (1 — x) ln(l — x)). 



(19) 



In the limit of low occupancy, which is of interest for us 
here, this expression is simplified as s{x) = x\ne/x. 

To make further progress, we assume that the interac- 
tion energy dependence on 6n can be approximated by a 
quadratic function. 



Eint{nA,nB) = -aSn^, 



N 



-SE{n), (20) 



where SE is the energy gain per adatom due to complete 
sublattice ordering, and n = (ua + nB)/2 is the average 
adatom concentration. Taylor expanding the entropic 
part of F in dn, we find 



y 



6 6 

{n -\- 6n) In h (f^ — 5n) In 



n + Sn 



n — Sn 



(21) 



Nn In- 



N 
2n~ 



Sn^ 



N 
12^ 



■Sn^ + OiSn^) (22) 



Using these results, we can write the effective free energy 
functional in the vicinity of the phase transition as 



AF = F{n- 
where 



Sn, n — Sn) — Fin, n) — N [a(5n^ 
SE{n) 



T 
2n 



T 
12n3' 



bSn^] , 
(23) 

(24) 



For the transition temperature, found from the condition 
a = 0, this gives 



Tc = 2SE{n) 



(25) 
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The values of SE taken from Fig. 2 inset (see the main ical temperature values are in the range of hundreds of 
text) are of the order of several hundred Kelvin. Thus, Kelvin, 
even for the adatom coverage as small as n « 0.03, crit- 



